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Abstract 

ON | All gasdynamical models for the evolution of gaseous content of galaxies assume 

■ that cooling from the hot, virialized phase to the cold phase occured through some 

sort of thermal instability. Subsequent formation of colder clouds embedded in 
the hot, rarefied medium is a well-known process appearing in many astrophysi- 
cal circumstances and environments. The characteristics of the condensed clouds 
00 . depend on the relevant timescales for cloud formation and disruption due to ei- 

CN ! ther collisions or one of the operating instabilities. In this paper, the importance 

of the Kelvin-Helmholtz instability is investigated for the clouds forming in huge 
gaseous haloes of L* galaxies. Recent treatment of this problem by Kamaya (1997) 
is extended and a more realistic cooling function employed. Results show that the 
Kelvin-Helmholtz instability proceeds effectively on the same timescale whether 
[~^« ■ we account for self-gravity or not. This has multiple significance, since these ob- 

, jects may have been seen as high-column density absorption line systems against 

the background QSOs, and probably represent the progenitors of the present-day 
globular clusters. 
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2 1 Introduction 

CO 

This is the first in a series of papers intended to discuss stability of metal-poor clouds in 
pressure equilibrium with the quasi-hydrostatic hot gaseous corona of (proto)galaxies. It 
seems that a sort of consensus exists that a state of virialized gaseous haloes located in 
the potential wells created mainly by the dark matter preceeded the formation of galactic 
subsystems as we know them today (Gott & Thuan 1976; White & Rees 1978; Rees 1978; 
Miyahata & Ikeuchi 1995; Mo & Miralda-Escude 1996). In the transition from this early 
stage to later steps necessary for the substructure formation, the major role certainly 
belonged to the process of cooling of gas. Various cooling scenarios were proposed, but 
it seems clear that the onset of thermal instability in a hot, virialized plasma led to 
creation of a two-phase medium, in which cold clouds are in pressure equilibrium with 
the ambiental, collisionally ionized medium. Cold clouds may be photoionized by an 
ionizing background. In the gravitational potential of the entire halo (i.e. baryonic + 
non-baryonic matter), cold clouds will certainly tend to fall toward the halo center. 

It is important to note that these physical processes are relevant not only for unob- 
servable early stages of galaxy evolution, but are operational at later epochs either, as 
in the Lyman-limit and metal absorption systems and, most probably, at least some of 
the lower column density Lya forest absorbers residing in extended haloes of luminous 
galaxies (e.g. Chen et al. 1998). Thus, a related motivation for undertaking research in 
this direction is the possibility to account for some fraction of the observed population 
of high-column density QSO absorption systems (Mo & Miralda-Escude 1996). 
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2 Field length in galactic-sized haloes 



The thickness of the boundary-layer where the phase transition of plasma is expected 
to be close to the so-called Field length, which is determined from the balance between 
the thermal conduction and radiative cooling (Field 1965; Aharonson, Regev & Shaviv 
1994). It can be written as 



A, 



V <A(T h ) ' 

where k is the " classical" coefficient of thermal conduction. The cooling function for the 
dominant free-free cooling is calculated in detail in Sutherland & Dopita (1993). 

The temperature dependence of the coefficient of thermal conduction can be, in the 
regime considered, approximated as (Cowie & McKee 1977; Kamaya 1997) 

k(T) = k T 2 - 5 , (2) 

where the constant is given as k = 5.6 x 10 -7 erg s -1 K _1 cm -1 . From Eqs. ([I]) and 
(H), explicitly taking into account the galactocentric distance of the considered cloud, we 
obtain 

7.5 x 10" 4 



cm. (3) 



n h (r) \| A[T h (r)] 

We see that the Field length is estimated from the physical condition (T h , n h ) of the 
hot plasma component. Its physical meaning is that the cold component will eventually 
evaporate (the timescale is discussed in detail in the Cowie & McKee [1977] paper) if its 
size is smaller than A^. On the other hand, the hot component will tend to condense on 
the cold cloud surface if its size is larger than A^. Now, our task is clearly formulated: 
to establish what is the value of A^ at various galactocentric distances in a particular 
gaseous halo model. 

The density and temperature profiles in the adiabatic model (Mo 1994; Mo & Miralda- 
Escude 1996) are given by 
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r \ 2 



and 



Ph (r) = p h (r c ) 1-iHn- , (4) 
\ r c / 

T h (r) = T h {r c ) (l - K\n -^) . (5) 
Notation is the following: K is a dimensionless constant equal to 

K 2 fxV[ 

bk B T h (r c y 

k-g, being the Boltzmann constant, fi average mass per particle, and Th(r c ) is temperature 
at the cooling radius, by assumption equal to the virial temperature: Th(r c ) = /zV c 2 /2k B 
(see Waxman & Miralda-Escude 1995 for an interesting discussion). For a typical V c ~ 
250 km s _1 corresponding to a L* galaxy, and a metallicity Z = 0.3 Z & this temperature 
is Th(r c ) = 4.89 x 10 6 K, and the constant K has the value K — |. The density at the 
cooling radius is obtained by requiring the cooling time at r c to be equal (at least) to t m : 

5pk B T h (r c ) 

Ph{r c ) = oX ,r r( yn , (7) 
2A[T h (r c ) \t m 
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Figure 1: Comparison between analytic approximation of Kamaya (1997), represented by 
dashed line, and more exact numerical calculations of Sutherland & Dopita (1993) for two 
different metallicities (solid lines) in the framework of adiabatic halo models. Evaluation 
is performed at fixed galactocentric distance of R = 100 kpc. 

where A[T^(r c )] is the cooling rate evaluated at the temperature of hot phase at the 
cooling radius. 

Kamaya (1997) has used a simple analytical approximation for the cooling function 
in the bremsstrahlung regime (see also Lepp et al. 1985) 

A(T) « 3.0 x 10" 27 T5 erg cm 3 s" 1 . (8) 

In Fig. 1, we see the difference between this form of the cooling function and the 
more precise numerical values of Sutherland and Dopita (1993). Cooling function for two 
relevant metallicities (0.1 Z & and 0.01 Z & ), in units of 10~ 23 erg cm 3 s _1 , is represented 
by solid lines, and analytic approximation of the eq. @ is shown as the dashed line. 
The comparison is made within framework of the adiabatic halo model (although the 
same applies to any other physical situation) at the fixed distance R = 100 kpc from the 
center of the halo characterized by the circular velocity V c \ r . The distance of 100 kpc is 
chosen because it is far enough from the region of the disk and non-stationary phenomena 
associated with the disk-halo connection. Simultaneously, it lies within the virial radius 
for almost all halo masses, so that the approximation of quasistationary cooling of gas 
can be legitimately applied. 

We see that for circular velocities characteristic for masses of galaxies (like the 
Milky Way), the analytic expression in (||) is a moderately good approximation for low- 
metallicity gas. It remains good within a factor of 2 up to the largest galactic masses 
(characterized by V cir ^300 km s -1 ), but tends to systematically underestimate the cooling 
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Figure 2: Field lengths for clouds at various galactocentric distances in the adiabatic 
haloes for V C \ T = 220 km s -1 (solid line) and V c \ v = 100 km s _1 (dashed line). 

rate for smaller galactic haloes. This is especially pronounced in case of dwarf galaxies 
(50 < Vrf r < 100 km s _1 ), where the disagreement is drastic. 

Now we use Eq. ([3]) in conjunction with Eqs. (|]) and ([5]). The resulting Field lengths 
are shown in Fig. 2 as functions of galactocentric distances for the two characteristic values 
of circular velocity. We see that Field lengths in general are quite large; clouds usually 
envisaged in adiabatic two-phase halo models (Mo & Miralda-Escude 1996; Chiba & Nath 
1997) will evaporate through thermal conduction, on the time scale t ev ap, as discussed by 
Cowie & McKee (1977). 

It is immediately clear that large values of the Field length are a consequence of 
extraordinarily small density of the hot plasma. It is easy to see that this value must 
always be low, irrespective of whether we are dealing with huge protogalaxies like those 
in Kamaya (1997) picture, or galactic-sized hot haloes of somewhat more modest size. 
If we perform a thought experiment and distribute the entire baryonic mass of an L* 
galaxy, M B ~ 5 x 10 11 M Q (e.g. Fields, Freese & Graff 1998) over the halo of radius 
R ~ 100 kpc, maximal density at the edge of such halo is ~ 5 x 10~ 3 cm -3 , and since the 
adiabatic profile is quite shallow, this can not be very different from our model values, 
implying Field lengths of tens of kpc. The existence of collapsed structures, like disk 
stars or MACHOs, will further increase this lengthscale. 

However, if the clouds are contiguous structures of the sizes ~ 100 kpc (and still 
surrounded by the hot medium) inferred by Dinshaw et al. (1994, 1995), they will be 
stable against evaporation. On the contrary, they will grow through condensations of hot 
ambiental plasma. 
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3 Application of the Richardson Criterion 



K-H instability which will tend to shred the cold component will be stabilized by self- 
gravity if the so-called Richardson criterion 

J > 0.25, (9) 

is satisfied. It is important to emphasize that other means of stabilization (through 
magnetic fields, for example), are irrelevant of the Richardson criterion, so the breakdown 
of this criterion does not automatically imply the onset of K-H instability. In the above 
inequality, the Richardson number J is defined as 

g dp/dr g p c - p h 
p c (dV/dr) 2 p c V? T 



The gravity of the cloud at position r is denoted by g, and the density of cold clouds 
by p c . Density and velocity gradients, dp/dr and dV/dr are evaluated at the boundary 
layer between the two phases, and corresponding vectors are oriented in such a way as 
to point from the cold to the hot component (e.g. dp/dr < 0). In the second part 
of this relation, we have used the approximation of Kamaya (1997), in assuming that 
dp = p h — Pc, dr = \ F . For gravitational acceleration, we may write the classical relation 

GM C 

9 = ^r, (ii) 

where mass and radius of the cold cloud are denoted by M c and R c respectively. Finally, 
we have approximated the relative velocity between the cloud and hot medium by the 
circular velocity V cir . There are several reasons why we should prefer this value to the 
approximation used by Kamaya (1997) which reduces to the virial velocity of gaseous 
subsystem only. The presence of large quantities of dark matter (95% by mass according 
to our fiducial adiabatic model) justifies using its virial velocity as the relevant velocity 
scale within entire halo. Besides, the motions of clouds are likely to be bound from above 
by terminal velocity (e.g. Benjamin & Danly 1997), which is in any case smaller than the 
circular velocity (see also the discussion in Mo & Miralda-Escude 1996). 

Using all these approximations, we can find Richardson number for two-phase medium 
of protogalactic or early galactic haloes in the form: 

GM C \ F 

J = 




M c \( R r \ V A; 



lO 6 M /\lkpc/ VlOkpc 



(12) 



We have used the fact that adiabatic models and other two-phase pictures predict pressure 
equilibrium of the two stable thermal phases with p h / p c ~ 1CT 2 . Of course, it is not 
necessary to achieve high precision in any of the terms in order to see that the resulting 
value is much smaller than the critical value of 0.25. 
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4 Discussion 



We have redone the calculations of Kamaya (1997) intended to demonstrate the indepen- 
dence of the Kelvin-Helmholtz instability of cloud self-gravity We have improved that 
work in two major respects: 

1. We have worked within the context of a specific, adiabatic halo model of Mo & 
Miralda-Escude (1996), instead of postulating appropriate temperatures and den- 
sities. 

2. The more realistic cooling function based on numerical calculations of Sutherland 
& Dopita (1993) was used instead of the analytic approximation. 

The main conclusion, however, remains the same: inclusion of self-gravity in the models 
of cold clouds formed by thermal instability does not change their subjection to the 
Kelvin-Helmholtz instability. This conclusion is valid not only for two-component model 
of protogalaxies of Miyahata & Ikeuchi (1995) and Kamaya (1996), but also for adiabatic 
models of galactic haloes of Mo (1994) and Mo & Miralda-Escude (1996) intended to 
explain the Lyman-limit absorption in QSO spectra. 

It is of foremost importance to model physical conditions in the early stages of galactic 
history, since that would enable us to understand the origin of the well-known features of 
galaxies observable today. To this end, it is necessary to achieve better contact of several, 
seemingly distinct, fields of astrophysical research. In forthcoming papers of this series, 
we shall discuss other physical processes affecting cold condensations in galactic-sized 
haloes and limiting their lifetimes, like the Jeans instability and collision timescales. 
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